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Abstract 

The atmospheric neutrino data collected by the Super-Kamiokande experi- 
ment span about four decades in neutrino energy E, and are thus appropriate 
to probe the energy dependence of the oscillation wavelength A associated to 
z^/i ^ flavor transitions, when these are assumed to explain the data. Such 
dependence takes the form oc i?" in a wide class of theoretical models, 
including "standard" oscillations due to neutrino mass and mixing (n = —1), 
energy-independent oscillations (n = 0), and violations of the equivalence 
principle or of Lorentz invariance (n = 1). We study first how the theoretical 
zenith distributions of sub-GeV, multi-GeV, and upward-going muon events 
change for different integer values of n. Then we perform a detailed analysis 
of the Super-Kamiokande data by treating the energy exponent n as a free 
parameter, with unconstrained scale factors for both the amplitude and the 
phase of <-> I'r oscillations. We find a best-fit range n = —0.9 it 0.4 at 
90% C.L., which confirms the standard scenario (n = —1) as the dominant 
oscillation mechanism, and strongly constrains possible concurrent exotic pro- 
cesses (n 7^ — 1). In particular, we work out the interesting case of leading 
standard oscillations plus subleading terms induced by violations of special 
or general relativity principles, and obtain extremely stringent upper bounds 
on the amplitude of such violations in the (i'^, Ur) sector. 

PACS number(s): 14.60.Pq, 96.40.Tv, ll.30.Cp, 04.80.Cc 
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I. INTRODUCTION 



The recent atmospheric neutrino data from the Super-Kamiokande (SK) experiment 
can be beautifully explained through flavor oscillations generated by nonzero neutrino 
mass and mixing j^J^ in the z/^ ^ Vr channel [^. Such interpretation is consistent with 
all the SK data, including sub-GeV e-like and /i-like events (SGe, /i) [Q, multi-GeV e-like 
and /i-like events (MGe,/i) 0, and upward-going muon events (UP/x) 0. A combined 
analysis of the 33 kTy SK data sample can be found in p[. The oscillation hypothesis 
has been strengthened by the latest (preliminary) 45 kTy SK data sample [pl,|lO|, and is 
also consistent with independent atmospheric neutrino results from the MACRO []TT| and 



Soudan-2 |]T2[ experiments, as well as with the finalized upward-going muon data from the 
pioneering Kamiokande experiment []T3|. 

Establishing ^ oscillations generated by nonzero v mass and mixing as the "stan- 
dard" interpretation requires, however, further data and analyses. Basically, the following 
three aspects should be clarified: (1) the periodicity; (2) the flavors; and (3) the dynamics. 

So far, the periodicity of the z/^ oscillation pattern has not been experimentally observed 
in the neutrino energy [E] or pathlength (L) domain, and it is unlikely to emerge from the 
SK lepton distributions, largely smeared in energy or angle. Although specific nonperiodic 
scenarios, such as neutrino decay [T^, can be indirectly excluded by careful analyses of SK 
data |[r5|jl6[| , the direct observation of a periodic disappearance pattern of z/^'s remains an 
important goal for future atmospheric and long-baseline neutrino experiments. 

Besides periodicity, one should also identify unambiguously the flavor(s) of the oscillating 
partner(s) of z/^'s because, in principle, all oscillation channels into active or sterile neutrinos 
might be open [y^ z/^, z/e, z/^ v^. While the amplitude of possible tran- 

sitions is bound to be small [H , one cannot exclude ^ oscillations with large amplitude 
with present SK data P,p!0|. However, based on the fact that different oscillation channels 
induce somewhat different energy-angle lepton distributions, there are good prospects for a 
significant (> 3(t) discrimination of z/^ Vr from z/^ Us with future SK data 



In this work we assume that both the periodicity (i.e., the existence of an oscillation 
length A) and the oscillating flavors (z/^, v^-) are established, and we rather focus on the 
third aspect to be clarified: The dynamical origin of z/^ ^ oscillations. The standard 
oscillation dynamics, involving a nontrivial 2x2 neutrino mass matrix, leads to a well-known 
energy dependence of A, 



A ^ oc £^ (standard) 



(1) 



However, possible nonstandard neutrino interactions or properties can also generate (or co- 

An incomplete list of possibilities include violations 
flavor- changing neutral currents (FCNC) (see 



exist with) z/^ ^ v-^ oscillations [|l8ip!9| 
of the equivalence principle (VEP) E 



2^ and refs. therein), neutrino couplings to space-time torsion fields p3[, neutrino inter- 



actions through charged scalar particles nonrelativistic heavy neutrinos violations 
of Lorentz invariance ( VLI) [0,^] and of CPT symmetry ||2^ . In several such models the 
energy dependence of A takes a power-law form p9| . 



A ""^ oc E"" (n 7^ — 1, nonstandard) 



(2) 
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Although models with exotic dynamics for z/^ ^ Vr oscillations do not survive Occam's 
razor, they might survive experimental tests! Effective tests must cover the widest possible 
energy range, as evident from Eqs. ([1|) and (|^). Concerning atmospheric neutrinos, pre- 
SK data analyses covered only about two decades in energy (i.e., the so-called contained 
events, E ~ 0.1-10 GeV), and were compatible with several nonstandard scenarios, in 
particular with the VEP hypothesis (corresponding ton = 1) ||30|j3l[| . An interesting post-SK 



analysis, covering a slightly more extended energy range [^, does not appear to discriminate 
significantly the three cases n = and n = ±1 examined. However, as observed in |jl5|,|16[, 
a much longer "lever arm" in the energy domain is provided by the inclusion of partially 
contained and upward-going muon events (up to ~ 10'^ GeV), thus providing a powerful 
tool to test exotic scenarios (which, indeed, appear to be disfavored in general [|T5[]). 



In this work we assess quantitatively the situation for z/^ ^ Uj- models with power-law 
energy dependence of the oscillation length (A~^ oc £'"), including, as relevant subcases, 
standard mass-mixing oscillations (n = —1) and violations of the equivalence principle or 
of Lorentz invariance (n = 1). We obtain two basic results: (1) The 90% CL. range of n is 
determined to be n = —0.9 ± 0.4, which is in striking agreement with standard oscillations, 
and excludes all n 7^ —1 models as dominant sources of Ur transitions; (2) assuming 

the n = 1 case as a subleading mechanism (coexisting with leading standard oscillations), 
we place stringent upper bounds on its amplitude. Such bounds can be interpreted as upper 
limits to violations of special or general relativity principles. 

The plan of the paper is as follows. In Sec. II we analyze models with generic oc E"^ 
behavior, and constrain n through fits to SK atmospheric u data. In Sec. Ill we consider a 
more complicated case with two coexisting sources for oscillations, namely, neutrino masses 
and violations of relativity principles. We summarize our results in Sec. IV. 



II. ANALYSIS OF MODELS WITH oc E"" 

In this section, we first review some oscillation models with A depending on E through 
a power law. Then, independently of specific models, we study the phenomenology of 
atmospheric z/'s, and constrain the energy exponent n through detailed fits to the SK data. 



A. Review of models 

Typical two-flavor hamiltonians for z/^ u^- oscillations predict a flavor transition prob- 
ability of the form 

Pii^^, ^ Vr) = sin^ 2^ sin^inX-^L) , (3) 

where ^ is the rotation (mixing) angle between the flavor basis and the basis where the 
hamiltonian is diagonal, L is the neutrino pathlength, and A is the neutrino oscillation 



length. In several cases of interest, A depends on the neutrino energy E as p9| , |32 | 

A"^ oc E" (4) 
with the exponent n taking integer (positive or negative) values. 
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For instance, the cases n = — 1, and 1 arise, respectively, in the presence of neutrino 
interactions mediated by scalar, vector, and tensor fields [^. Specific models include: 



A 



-1 



At^E 

TT 

E6v 



27r 
5h_ 



n 



n 



n 



n 



n 



—1, standard oscillations [2]) , 
1, violations of equivalence principle 
1, violations of Lorentz invariance ||26|| ) 



0, violations of CPT symmetry |28 



0, nonuniversal coupling to a torsion field [2^) 



(5a) 
(5b) 
(5c) 
(5d) 
(5e) 



Standard oscillations [Eq. (|5aD ] are simply generated by nonzero neutrino masses (Am^ = 
m\ — ml), with a mixing angle ^ [Eq. (^] usually denoted as 6^23, ip, or simply 9 (we adopt 
the symbol ip in the following, as we use 6 for the lepton zenith angle). In this case, the 
energy exponent is = — 1. 

Equation (|5b| ) refers to possible violations of the equivalence principle (VEP) pO| , 
namely, to nonuniversal couplings of neutrinos to the gravitational potential 0. The differ- 
ence in such couplings is usually denoted by A7, and the mixing angle by 9g- The potential 
seems to be dominated by the local supercluster (0 ~ 3 x 10~^ [0); but ambiguities in its 
definition suggest to use the product 0A7 as a single, dimensionless free parameter, rather 
than (p and A7 separately (see PT|P^ , P3[] for recent discussions). The energy exponent of 
VEP-induced oscillations is n = 

Equation (^ refers to possible violations of Lorentz invariance (VLI), namely, to asymp- 
totic neutrino velocities different from c The parameter 6v represents the speed differ- 
ence in units of c. The mixing angle is usually denoted as 6^. Since the energy exponent n 
is the same (= 1) in both the VEP and the VLI scenario, such mechanisms are phenomeno- 
logically equivalent through the substitutions I0IA7 ^ 6v/2 and 9g 0^ p7| . 

Equation (^) refers to possible violations of the CPT symmetry through a more general 
class of Lorentz- violating perturbations p8[, 5h being proportional to the CPT-odd hamilto- 
nian producing energy independent (n = 0) oscillations. It is interesting to notice that the 
most general Hamiltonian considered in |^ encompasses the three scenarios with n = — 1, 
and 1. 

Finally, Eq. (Q refers to possible nonuniversal couplings (Afc 7^ 0) of neutrinos to a 
space-time torsion field of strength Q which also produce energy-independent oscilla- 
tions. 

In principle, several oscillation mechanisms with different n's might occur at the same 
time, with corresponding complications in the analysis. In the next two subsections, however. 



^An alternative, string-inspired VEP mechanism leading to an energy exponent n = — 1 rather 
than n = 1, has been recently considered in |37|. Its phenomenology would be indistinguishable 
from the standard far as ^ u^- oscillations are involved. We do not consider such VEP 

scenario in this work. 
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we consider only one mechanism at a time. We will discuss the interesting case of coexisting 
n = —1 plus n = +1 oscillations in Sec. III. 

A remark is in order. Not all "exotic" models for ^ transitions can be parametrized 
as in Eqs. (H) and (|). An important exception is represented by FCNC- induced oscillations 
(see and references therein). In fact, although the FCNC oscillation phase is energy- 
independent (n = 0), it is proportional to the column density of electrons rather than 
to L. Therefore, FCNC-induced oscillations deserve a separate analysis |22|,|T3| and will be 
considered in a future work. However, some features of the n = case also apply qualitatively 
to the FCNC case. 

Finally, we note in passing that the analysis of nonstandard energy dependences for u 
oscillations has some correspondence in the neutral kaon system, where the role of A is played 
by the effective K°-K° mass difference (see and references therein). 



B. Model-independent analysis 

In this section we do not commit ourselves to any specific model, and rather analyze the 
phenomenology of the oc dependence in the most general way. We assume that the 
^ oscillation probability takes the form 

P(.„«..) = „.=in^(,3 — — ) , (6) 

where a is an overall scale factor (0 < a < 1) for the oscillation amplitude, f3 is an un- 
constrained scale factor for the oscillation phase, and n is a free exponent (not necessarily 
equal to —1, 0, or 1). The units in Eq. (|^) have been conveniently chosen on the basis of 
the current Super-Kamiokande data, which suggest an oscillation length of 0(10^ km) for 
neutrino energies of 0(1 GeV). The standard oscillation case is recovered by taking n = — 1, 
a = sm'^2^p, and /3 = 1.27AmV(10-3 eV^). 

Equation (^) is used to calculate the observable rates of sub-GeV, multi-GeV, and upward 
through-going muon events in SK as a function of the lepton zenith angle 6, using the same 
detailed and accurate approach as in Refs. f^JT^, to which we refer the reader for technical 
details. Here we just remind that the parent neutrino distributions for sub-GeV (SG), 
multi-GeV fully contained (MG FC), multi-GeV partially contained (MG PC), and upward 
through-going (UP) muon events are roughly distributed in the ranges 0.1-3 GeV, 1-10 
GeV, 1-10^ GeV, and 10-10^ GeV, respectively, thus covering about four decades in E. 

Figure 1 shows the expected zenith distributions of SG/i, MG/i (FC+PC), and UP// 
events in Super-Kamiokande for maximal mixing (a = 1) and for three representative val- 
ues of P: dotted line, /5 = 0.2; dashed line, /3 = 1; solid line, /3 = 5 (on color printers 
such hues appear in green, red, and blue, respectively). In each bin, the muon rates fi are 
normalized to the expectations /xq in the absence of oscillations, so that deviations from 
the no-oscillation case (/i//io = 1) are immediately recognizable. Since we are considering 
pure z/^ ^ Ur oscillations (e/co = 1 always), the electron rates are not shown. The effect 
of taking a < 1 can be approximately described by reducing proportionally the amount of 
muon disappearance in each subfigure. The theoretical expectations are affected by rela- 
tively large and strongly correlated uncertainties (not shown), mainly related to the overall 
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normalization of the atmospheric neutrino flux The 45 kTy SK data (dots with 

statistical error bars) are superposed to guide the eye, although no fit to the data is implied 
by this figure. 

In Fig. 1, from top to bottom, the energy exponent n takes the values n = —2, —1, 0, 
and 1, corresponding to a dependence of the kind L/E'^, L/E, L, and L-E ioi the oscillation 
phase. Such four cases are described in the following. 

Case L/E"^. This case {n = —2) does not correspond to any known model, and is used 
only to extend the study to a power-law energy dependence faster than in the standard case 
(n = —1). The muons appear to be significantly suppressed at the lowest (SG) energies, 
the more the larger (3. An excessive suppression of SG muons is avoided only by keeping 
P ^ 1. On the other hand, such values of P are too low to produce a significant suppression 
of MG muons, which rather prefer /5 ^ 5. Therefore, one expects a "compromise" between 
underestimated SG/x rate and overestimated MG/i rates for P in the range ~ 1-5. UP/^ 
events are basically unsuppressed, due to their high energy. Reducing a would be of no help 
in reducing the conflict between SG and MG muon expectations. 

Case L/E. This is the "standard" oscillation case, which, as well known, provides and 
excellent description of the SK data for Am^ ~ few x 10~^ eV^ Therefore, it is not 

surprising to see in Fig. 1 that, for values of (3 in the range ~ 1-5, all the muon data are 
well reproduced. Higher values of P would suppress SG muons too much, while lower values 
would not produce enough suppression at higher energies (MG/x and UP/i samples). 

Case L. In this case, the z/^ survival probability is energy independent, and the differences 
in the muon suppression patterns among the SG, MG and UP muon samples are mainly 
due to the different angular smearing. At low (SG) energies the smearing is very effective 
and there is little difference between the three curves, while higher-energy, MG muons are 
more discriminating. Values of P around unity (dashed curve) seem to be preferred by 
SG+MG data, and produce a significant (~ 1/2) suppression below the horizon. UP muons 
have a different suppression pattern, which is highly correlated with the parent neutrino 
direction, and follows closely the variations of P^^ with L (dashed and dotted curves), until 
the oscillations are so fast to be unresolved within the bin width (solid curve). The expected 
UP/i suppression appears to be larger than suggested by the data, unless a is taken to be 
nonmaximal; however, for a < 1 the relatively good description of SG+MG data would be 
spoiled (not shown). 

Case L-E. In this case, the z/^ survival probability rapidly approaches the average value 
1/2 as the energy increases. Therefore, although SG+MG data can be described relatively 
well with a = 1 and P ~ 0.2, the expected UP/z rates are too low in any case. As in the 
previous n = scenario, taking a < 1 would help to reproduce the UP/i data, but would 
worsen the description of SG+MG data. No satisfactory compromise can be reached. 

In all the four cases, it can be seen that SG (or even SG+MG) data alone do not 
discriminate strongly among the various scenarios, since they do not probe the full energy 
range explorable by SK. This might explain why the SG+MG analysis in could not 
significantly distinguish the three cases n = —1, 0, and 1. The inclusion of UP/i data 
extends the range up to ~ 10^ GeV, and provides a very important tool to probe the 
energy dependence of A |T^. 

In conclusion, from the examination and the comparison of the oscillated muon distribu- 
tions at different values of n shown in Fig. 1, the case n = —1 emerges as a good description 
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of the SK data at all energies, while for n ^ —\ the patterns of muon suppression at low 
(SG), intermediate (MG) and high (UP) energies appear to be in conflict with the data. 



C. Fits to the Super-Kamiokande data 

The qualitative understanding of the muon distributions at different values of n (previous 
subsection) can be improved by performing quantitative fits to the 45 kTy preliminary SK 
data P,p!0[|. We use a approach that, as described in P], takes into account several sources 
of correlation among the systematics affecting the theoretical predictions. Even if i/g's do 
not participate to oscillation in the scenarios considered here, we have included the SG and 
MG e-like data in the analysis (for a total of 30 data points), since they play an important 
role in constraining the overall normalization uncertainty (see also P,p!6|). 

Figure 2 shows the best-fit zenith distributions of muons for the four cases considered in 
Fig. \ [n = —2, —1, 0, and 1), as obtained by leaving a and /3 unconstrained. The values 
of a, /3, and at the best-fit points are reported in the top part of the figure. The LjE 
(n = —1) case provides an excellent fit to the data (XminZ-^DF = 20.3/28), while all the other 
cases do not provide a good description of at least one data sample (SG/i, MG/i, or UP/i). 
In particular, for n = —2 there is insufficient up-down asymmetry of MG/i's and no slope of 
UPyLi's; for n = none of the zenith distributions is correctly reproduced; for n = 1 there is 
a too strong and fiat suppression of UP/x's. 

In Fig. 3 we present the curve as obtained by taking also as a free parameter, besides 
a and (3. Although noninteger values of n may not be related to any realistic oscillation 
dynamics, this exercise is useful to see how accurately n is determined through the SK data. 
The result is striking: 

n = -0.9 ± 0.4 at 90% C.L. (7) 

(corresponding to ~ xLm — 6-25 for A^df = 3). This narrow range for n is perfectly 
consistent with standard {n = —1) neutrino oscillations, and inconsistent with any other 
integer value of n. 

Given the importance of standard z/^ ^ Ur oscillations, we show in Fig. 4 the updated 
limits on the oscillation parameters Am^ and sin^ 2ip. We find the best fit at Am^ = 
2.8 X 10~^ eV^ and maximal mixing. The bounds in Fig. 4 are in good agreement with the 
latest full data analysis from the SK collaboration | |10| ]. 

In conclusion, standard oscillations (A~^ oc E~^) are strongly favored as the dominant 
mechanism for the z/^ ^ flavor transitions of atmospheric neutrinos in SK. Alternative 
mechanisms of the kind A^^ oc -E" (with n ^ —1) cannot be the dominant source of muon 
disappearance in SK. In particular, violations of special or general relativity principles (VLI 
or VPE, leading to oc E) cannot explain the bulk of SK atmospheric u data. Therefore, 
if n 7^ — 1 oscillations occur in nature, they can only be subleading processes with small 
amplitude, coexisting with leading, large-amplitude n = —1 standard oscillations. Such 
results generalize and refine previous indications that SK data could disfavor some exotic 
models [ITH . 
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III. CONSTRAINTS ON VIOLATIONS OF RELATIVITY PRINCIPLES 



In this section we consider a more complicated case, characterized by leading n = —1 
oscillations plus subleading n = +1 oscillations, possibly generated by violations of relativity 
principles (YEP or VLI). We show that the fit to the SK data is not improved with respect 
to the case of standard z/^ Ur transitions. As a consequence, we derive upper bounds on 
such violations. A brief review of the theoretical formalism precedes the phenomenological 
analysis. 



A. Violations of relativity principles: Formalisms 



The theory and phenomenology of neutrino violations of the equivalence principle (YEP) 
20| have been investigated in a number of papers, including studies of the solar u deficit 



30| , |3llj3^ , |35|j39| - ^] , of the atmospheric u anomaly p9|-|3^, of oscillation searches at short 
baseline and long baseline [^,^,0 accelerator facilities, and of double beta decay 

47 1 . Given the phenomenological equivalence of violations of Lorentz invariance (YLI) 
and of the equivalence principle (YEP) p^, neutrino oscillation searches can be generally 



interpreted as tests of fundamental principles of both special and general relativity, with a 
sensitivity at levels below 10~^° (see, e.g., P5|PT| ). Here we focus on the case of YEP or YLI 
induced oscillations coexisting with standard oscillations. Mixed scenarios of this kind have 

mm. 



been considered, e.g., in 



but, to our knowledge, they have not been discussed 



on the basis of Super-Kamiokande atmospheric u observations so far. 



In the presence of several concurrent processes leading to 



oscillations, the global 



Hamiltonian H is the sum of several 2x2 matrices Hn, which can be diagonalized through 
separate rotations (with angles of the flavor basis (z/^, Ut). As far as models of the kind 
oc are concerned, the rotation angle ^, which diagonalizes the total hamiltonian, is 



related to the oscillation length A through equations of the form 28 



TrA^Sin 2^ 



^c„ sin2a^"e''"' 



vrA"^ cos 2^ = Y,Cn cos 2^™ 



(8a) 

(8b) 



where the coefficients c„ parametrize the strength of each oscillation mechanism. In general, 
only one of the complex phase factors e^^" can be rotated away, the others being physically 



observable P8| , |3T| . For any given choice of the parameters (c„, //„), one has to derive the 
values of A and ^ from the previous equations, and insert them in Eq. (^ to get the flavor 
transition probability. 

In the specific case of standard+YEP {n = —l(Bn = +1) oscillations, Eqs. (|8aD and (|8bD 



can be rewritten as [31 



ttA"^ sin 2^ 



1.27 sin 2^ + 5.07 E sin 2^^ e'"^ 



E 



ttA"^ cos 2^ = 1.27—— cos2^/' + 5.07 



E 



10-21 

I0IA7 
10-21 



E cos 200 , 



(9a) 
(9b) 
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where the following units have been used: [Am^] = 10~^ eV^, [L] = [A] = 10^ km, and 
[E] = GeV. The same equations formally apply to violations of Lorentz invariance, modulo 
the replacements I0IA7 6v/2 and 9g 0^ Notice that the oscillation phase n\~^L, 
to be inserted in Eq. (|^), is proportional to the geometric average of the right hand sides 
of the above equations, so it will contain, besides the standard term (oc E^-^) and the VEP 
term (oc E), also an energy- independent interference term. The oscillation amplitude sin^ 2^ 
also acquires a nontrivial dependence on the neutrino energy through Eq. (|^). 



B. Constraints from Super-Kamiokande data 



Before performing a detailed fit to the SK data in the standard+VEP oscillation scenario, 
let us derive some qualitative bounds on the magnitude of possible VEP terms. According 
to the conclusions of Sec. II, we expect that the second (VEP) term on the right hand of 
Eqs. ( pa| ) and (pB[ ) should be typically much smaller than the first (standard) term, namely 

^•^^ 10-21 GeV 10-3 eV^ E ' ^^"^ 

or equivalently (for Am^ ~ few x 10^'^ eV^): 

|0|A7«lO--(^)' . (11) 

Therefore, the constraints to VEP effects should be stronger in the highest-energy SK data 
sample (UP/x events). Since the parent neutrino energy spectrum for UP/i's is peaked around 
10^ GeV, the sensitivity to VEP-induced oscillations is expected to reach levels of 0(10-^^) 
in the parameter |0|A7 (or, equivalently, in the parameter Sv/2 for the VLI case). Of course, 
such sensitivity depends somewhat on 9g- 

In order to get some insight about the 6'G-dependence of the expected constraints, let us 
consider the extreme values 9g = and 6g = 1^, and take (Am^, sin^ 2?/') at their best-fit 
values (2.8 x 10"^ eV^, 1). We also fix e*'' = 1 for simplicity. Then P^^ takes a simple form, 

-^sin2(Av/rT^) , ^G = 0, 

^M. = <! (12) 

{A{1+x)) ,^G = vr/4, 



sin^ 



where A = 3.56 L/E, x = 1.42 E^j^l A7/IO-21, and the units are [E] = GeV and [L] = 
10^ km. The standard case {P^^^ = sin^ A) is recovered for a; = 0. In the UP/x event sample, 
where the VEP effect is larger, the value of A is typically small (^ 0.6), and one can easily 
check numerically [from Eq. (0)] that the difference P^^ — P^**^ grows more rapidly with x 
for 6g = 7^/4: than for 6g = 0. Therefore, we expect a higher sensitivity to x (i.e., to I0IA7) 
at larger 6g- Moreover, at small values of both A and x it turns out that P^r — < 
(> 0) for Og = {6g = 7i"/4), so that the muon rates should be less (more) suppressed than 
in the standard oscillation case. 

Figure 5 illustrates quantitatively the previous considerations, by showing the VEP effect 
(added to best-fit standard oscillations) on the SK muon distributions for two representative 
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cases: (i) Oq = and |0|A7 = 1.5 x 10"^"^ (solid line, blue in color); (ii) 9g = 11/4: and 
I0IA7 = 2.0 X 10~^^ (dotted line, green in color). We anticipate that the such values are 
close to the border of the parameter region excluded by SK. The standard oscillation curves 
(I0IA7 = 0) are also shown for reference (dashed lines, red in color); they are identical 
to the best fit curves for standard oscillations (n = —1) in Fig. 2. As expected from the 
preceding discussion, the VEP effect is manifest at high energies (UP/x sample), and the 
expected deficit of UP/x's is more pronounced for 6*0 = 7r/4 than for 9g = 0. 

The next step is to perform a analysis of the standard+VEP scenario. We take for 
the moment (Am^, sin^ 26*) = (2.8 x 10~^ eV^, 1) and e*'' = 1, while leaving the parameters 
(|0|A7,^G') free. We find an important result, that strengthens the conclusions of Sec. II: 
The fit is never improved in the presence of VEP effects, as compared with the value 
Xmin = 20.3 derived in Fig. 2 for the standard case. Thus, not only VEP-induced oscillations 
cannot be the leading mechanism underlying the SK observations, but also there is no 
indication in favor of subleading VEP oscillation terms. As a consequence, we can place 
well-defined upper bounds on violations of relativity principles in the (z/^, z/^) sector. 

Figure 6 shows the 90% and 99% CL. limits on I0IA7 as a function of the VEP mixing 
parameter sin^26'G- The same limits apply to the neutrino asymptotic speed difference 
Sv/2, as a function of the VLI mixing parameter sin^26'„. As expected from the discussion 
at the beginning of this subsection, the limits obtained in Fig. 6 are roughly of O(10~^^), 
and become stronger as 6g increases. The stringent bounds in Fig. 6, together with the 
results shown in Fig. 3, represent our main contribution to the current understanding of 
atmospheric neutrino oscillations induced by violations of relativity principles in the (z/^, Ur) 
sector. Notice that such bounds pre-empt the region of VEP (or VLI) parameters explorable 
with proposed long-baseline accelerator neutrino facilities |]3l||.p| 

Finally, we have investigated the robustness of the bounds shown in Fig. 6 under vari- 
ations of the standard mass-mixing parameters. We have repeated the fit by varying Am^ 
and sin^ 2ip within the 90% CL. limits shown in Fig. 4, and also by taking negative values 
for A7, as well as generic values for the complex phase rj. We have found that, in any 
case, the value of xLin is not smaller than in the best-fit standard case. Therefore, stan- 
dard+VEP (or standard+VLI) oscillations never represent a better description of the SK 
data, as compared to best-fit standard oscillations. Under the above variations, the upper 
bounds shown in Fig. 6 are somewhat modified within factors of a few, but do not change 
qualitatively: they always become stronger as sin^ 29g increases. The most conservative 
upper bound (including negative A7 cases) turns out to be 

I0A7I < 3 X 10"^^ at 90% CL. , (13) 

independently of 9g- In the specific case 9g = 7r/4, the above bound can be lowered at least 
to ^ 10~^^. Analogous limits apply to the VLI parameters |5f |/2 and 9y. In particular, 

\6v\ < 6 X 10-24 at 90% CL. (14) 



^Notice that the values of L, E, and L • E probed in the UP^ sample by Super-Kamiokande are 
higher than in proposed long-baseline neutrino beams. 
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To our knowledge, the above limits to violations of special or general relativity principles 
are the strongest ever placed in the (z/^, u^) sector. They are valid under an assumption 
which is supported by the present data and appears likely to be corroborated in the future, 
namely, that standard i/^ 1/^ oscillations generated by ly mass and mixing represent the 
dominant mechanism underlying the Super-Kamiokande observations. 

IV. SUMMARY AND CONCLUSIONS 

Among the z/^ ^ z/^ models with oscillation length following a power-law energy de- 
pendence, standard neutrino oscillations generated by u mass and mixing are unique in 
providing a good description of the Super-Kamiokande atmospheric neutrino data, and are 
strongly favored as leading mechanism for (z/^, 1/^.) flavor transitions. Additional, subleading 

^ oscillations generated by possible violations of special or general relativity in the 
neutrino sector do not improve the agreement with the data, and must thus have a relatively 
small (or zero) amplitude. In particular, the fractional difference of asymptotic u velocities 
|5f 1/2 or of z/ couplings to gravity {(pAj] cannot exceed the value ~ 3 x 10"^'' at 90% C.L. 
for unconstrained neutrino mixing. The broadness of the neutrino energy range probed by 
Super-Kamiokande is crucial to obtain such strong limits. 
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Super-Kam zenith distributions of muons 
for oscillations with max. amplitude, a = 1 

P^^ = a - s\u(J * L/1 O'km * E7I GeV") 
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FIG. 1. Expected zenith distributions of SG^u, MG/x, and UP/i events in Super-Kamiokande 
for maximal mixing (a = 1) and for representative values of (3: dotted line, (3 = 0.2; dashed line, 
/3 = 1; solid line, /3 = 5 (on color printers such lines appear in green, red, and blue, respectively). 
In each bin, the muon rates ^ are normalized to the expectations /io in the absence of oscillations. 
The 45 kTy SK data (dots with statistical error bars) are superposed to guide the eye. From top 
to bottom, the energy exponent n takes the values n = —2, —1, 0, and 1. 
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FIG. 2. Best-fits to the zenith distributions of muons in SK for the four cases considered in 
Fig. l(n = — 2, — 1, 0, and 1), as obtained through a analysis of all the data (//-like and e-like 
events) with unconstrained a and /3. The values of a, /3, and at the best-fit points are reported 
in the top part of the figure. 
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FIG. 3. function from the fit to the SK data, assuming continuous vakies of the energy 
exponent n and unconstrained scale factors for the oscillation amplitude a and phase /3. The 
standard case (n = — 1) is very close to the best fit. 
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standard i^^^^Vr oscillations 
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FIG. 4. Updated bounds on the neutrino mass-mixing parameters for standard Vr os- 

cillations, as derived by our global analysis of all the Super-Kamiokande atmospheric neutrino 
data. 
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standard i^^.— >i^t oscillations at best fit 
+ violations of the equivalence principle 
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FIG. 5. Effect of subleading oscillations induced by violations of the equivalence principle (VEP) 
on the muon distributions in Super-Kamiokande. Standard mass-mixing oscillation parameters are 
taken at their best fit, and the complex phase e*'' is taken equal to 1. Dashed line (red in color): 
no VEP. Sohd line (blue in color): VEP with |(/>|A7 = 1.5 x 10~^^ and Qq = 0. Dotted hue (green 
in color): VEP with |(/)|A7 = 2 x lO^^e and = vr/4. 
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standard z^^^i^r oscillations at best fit 
+ violation of relativity principle(s) 



Violation of 
Equivalence 
Principle 

\(p\^y 



10 



-24 



10 



-25 



10 



-26 



SK 45 kTy: SGc,/i + MGe,^ + UPjj, (30 data) 

90% C.L. (Ndp = 2) 

99% C.L. 



Excluded 



- 10 



Violation of 

Lorentz 
Invorionce 

2 



-24 



- 10 



-25 



-26 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



10 



G, V 



FIG. 6. Bounds on the parameters that characterize violations of special or general relativity 
principles, assumed to generate subleading <-> v^- flavor transitions concurrent with standard 
(leading) neutrino oscillations. Best-fit neutrino mass- mixing values are assumed. See the text for 
further details. 
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